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We show explicitly that the one- and two-dimensional Hubbard model does not show phase sep- 
aration at any filling away from half-filling at T = 0. Apart from a single plausible assumption, 
only known exact results and symmetry properties for the one-band Hubbard model are used. Im- 

. . . plications for usage of the simple one-band Hubbard model for analysing the physics of cuprate 

superconductors are discussed. We also discuss models where electronic phase separation might 

O : generically occur, in the specific context of colossal magneto-resistance manganites. 

■ PACS numbers: 

o : 

■ The normal state properties of oxide superconductors are as perplexing as their high transition temperatures. 
^ ' Considerable experimental evidence has revealed that many high Tc superconductors have a regime way from half- 

filling in which electronic phase separation occurs. Such a phase occurs between an oxygen-poor and an oxygen-rich 
[ phase that goes superconducting. It is believed in some works, therefore, that the presence of electronic phase 
|~~] ^- separation (EPS) is connected with the physics of these doped materials. 

O [ Au et al and Su [3,4] have obtained rigorous results on phase separation in the 2D Hubbard model on a bipartite 
^ lattice. However, in [4], the question concerning the possibility of phase separation cannot be resolved at T=0. If 
EPS would indeed exist in the 2D Hubbard model at T = 0, one would suppose the precursor effects associated with 
this T = instability to manifest themselves in the physical response of cuprates, especially near the carrier doping 
driven T = Mott transition. For example, a singular compressibility is implied near such a quantum transition; this 
would necessarily affect the charge dynamics in a drastic way near the IMT. 

In this paper, we make use of an inequality due to Pitaevskii et al [5] to address this issue. Unlike in [4], which 
^ makes use of rigorous results on the 2D Hubbard model at T > for a square lattice, this inequality is valid for T=:0, 
Q ^ and contains information about the quantum fluctuations in the system. 

O . We start with a one-band Hubbard model defined on a square lattice with M sites in an external space-dependent 
\ magnetic field 

with a staggered external Zeeman magnetic field, —{h/2) crnio-e^'*''^', where the symbols have their usual mean- 
ings. We will consider eqn (1) with bot h signs of U. Inspite of intensive study, very few exact results are known for 
this model. It has been shown that there is no spontaneous magnetic order in one- and two dimensions at any finite 
(yQ ' temperature, as is expected from the const raints imposed by the Mermin Wagner theorem [6]. The exact solution of 
. the ID Hubbard model has been worked out by Lieb and Wu [7] using the Bethc ansatz. 

Pitaevskii's inequality, which provides direct information on quantum fluctuations of the system, reads [5] 
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{[AlAU){[B\BU)>\{[A\B])\^ (2) 

for any two operators A and B satisfying (0|A|0) = = (0|i3|0) (|0) is the ground state of H). Here, denotes 
the anticommutator of A and A^ and g is fixed by the condition 

gig-Ri = 1 if R_j connects sites on the same sublattice, and equals -1 otherwise. Eqn (2) then gives 

S^iq + g)S^{cD>[sir (3) 

where we used [S'+,5'-] = 2s^ and introduced the static structure factor S^{k) = {l/2N){[Sl{'k), S+{'k)] + ) along 
with the staggered order parameter 

4 = (l/iV)(E,s^e's-^') (4) 
for the antiferromagnet. Following [5], we consider the Cauchy-Schwarz inequality. 



^^(q) = J S^{q,uj)dLU <J J ijS^i(i,Lu)duj J ^^^^^dLJ (5) 
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and introduce an effective speed of spin excitations by 



2 .so _ Ju;S^{(l,u;)du; 



so that eqn (3) becomes 



si 



5^(q + g)> 7\^, (7) 

9Ce//(q)x-^(q) 



where X"'"(q) = / ^—^^^du is the transverse static response function. 
To proceed, we consider 

j uS^{ci,iv)cL, = (l/2iV)([5|(q), [i?,5+(q)]]> (8) 

with H given by eqn (1) on a square lattice. Explicit evaluation of the double commutator yields the result 

/ u;S-^{ci,u;)dw = ^(£k-q - £k)(nkT - «ki) + 2h{S^{-g)) (9) 
k 

where £k = Tiite^^'^'. Also , it is easy to show [6] that 

t ^(1 _ e^q-R.) = t ^(1 - cos(q • RO) (10) 

i i 

Using this in eqn (9) gives 

{[Slid), [H,S+{ci)]]) = Y^[q'R^t{ROP{Ri) + 2/i5o.(-g)] (11) 

i 

where -P(Ri) = (l/2)Ske~*'''^'(nk| — n^i). The i(Ri) are matrix elements of the one-electron operators between 
Wannier states, which fall off rapidly with Ri. So Si(l/2)J??i(Ri)P(Ri) = L is well defined. Also, Soz{-g,) = 
(l/A'')(S';j:(— g)). With these substitutions, eqn (11) becomes 

(l/2iV)([5t.(q), [H, 5+(q)]]) = il/2)[Lq' + 2hSo.{-g)] (12) 
Substituting into eqn (5) yields 



S^id) = I S^{ci,u)d^ < sj [^mW + 2/i5o.(-g))] (13) 

Using the definition of (g) , we have [5] 

Sq5^(q + g) < (l/2)(ArV2Ar) = N/A (14) 

Hence, 

(iV/4) >slY, [(i/2)^±(q)(L52 + 2/i5o.(-g))]V2 



(15) 



or 

H -1 

[5o.(-g)r<(Ar/4)'" 



E 



(x^/2)[V + 2/i5o.(-g)]V2 



(16) 



q 

Here, we used the fact that the perpendicular susceptibility X^{^) approaches a constant value for small (j[5]. Replacing 



the sum over q by an integral Sq ^ J (i-n)^ -i letting be the distance of the nea rest Bragg reflection plane from 
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the origin in q-spacc, we evaluate the integral to get 7(1) = In [(l/a)^/^ + (1 + in £> = 1, while in 7) = 2, 

we get 7(2) = [go(l + a)^/^ - a^/^] . Here a = ^^^^°|iy^" , so that for small h, we get 



[^o.(-g)F<(^)-0 ash^Oind=l (17) 



[So.i-S)? < (^) > as h ^ in d=2 (18) 

Hence, our analysis shows explicitly the absence of magnetic order in one dimension, but yields yields only an upper 
bound at T = in 2d. It is believed in the literature that the 2d Hubbard model at T = exhibits antifcrromagnctic 
long range order [7(a,b),8]. Together with the rigorous result on the spin-spin correlation function of the Hubbard 
model on bipartite lattices [11], our result uses only exact results and symmetries of the one-band Hubbard Hamiltonian 
on bipartite lattices to show that the existence of antifcrromagnctic LRO is not excluded in 2 dimensions. Wc 
emphasize that the derivation above does not rely on the applicability of the semiclassical spin wave approximation, 
but treats the dominant T = quantum fluctuations. We also notice that 5oz(— g, h) is a smooth function of h for 
small h. Wc shall make use of this fact in what follows. 

We now use the unitary particle-hole transformation for a bipartite lattice [9], Ci| — > Ci|, Cjj, — > e{i)Cl^ with 
e{i) = 1, i in A-sublattice, and equals —1 otherwise. The Hubbard Hamiltonian is transformed to 

H = -t Y,{ClCj, + h.c) - [/ ^ - {{h - U)/2) ^ n,, (19) 

ija i ia 

and 



He,t = -(m - {U/2)) J2{ni^ - na) (20) 

i 

where i7ea:t is the part of the transformed hamiltonian that looks like a Zeeman term describing the linear coupling 
of the fermions to a magnetic field /i = (/x — {U/2)). 

In the case of a singlet ground state, 2{S^S~) = {SfS^) = e{i)e{j)Cij w ith C,j > 0[11]. Under these conditions, 
\Sf \ < \S^{—g,h)\. Using this inequality, the eqn for |[S'o2(— g, /i)]| takes the form 

1 - C[lnl3]-^/^ < Pin') < 1 + C[ln(3]-^/^ (21) 

where j3 = 2 {qlLl2{2ij! - U')) in d = 1. Also, 

1 - C"[(2//3) + Cr"^ < p(m') < 1 + C"[(2//3) + C]-^/' (22) 

where ^ > 0, in (i = 2. Here C and C are constants, a nd the eqs (23) and (24) are valid for small {2iJ.' — U'). Also, 
i<T{f^i<T)i the elect ronic density per site. The above equations show that in one and two dimensions, for 
sm,all dopings. i.c for fi' close to U'/2, the density, p{fi), is a continuous fimction of p. According to the criterion fo 
r phase separation requiring a discontinuous variation in density with chemical potential [10], we conclude that the 
one-band Hubbard model near half-filling does not exhibit the phenomenon of phase separation in either one or two 
dimensions at T = 0. At exactly half- filling, the ground state of the Hubbard model is expected to be insulating both 
in one- as well as in two dimensions for any U/t on a bipartite 1 attice. This holds for our case of a two-dimensional 
square lattice, and is also the case for a one dimensional linear chain. At an arbitrary density away from the half-filled 
case, we define pseudospin operators 

J+ = eii)C}^Cl, Jr = [J+]t, Jf = (1 - - n,i)/2 (23) 

satisfying the usual pseudospin SU(2) algebra. Without loss of generality, we can set the magnetic field h to zero in 
what follows. Choosing A = (q + g) , and B = J~^{q), wc can show, using the Pitaevskii and Stringari inequality, 
that [J^] is a smooth, continuous function of {p — {U/2)) in one- and two dimensions at T = 0. We also reach the 
same conclusion by using an argument due to Shen et.al [11], based on the unitary mapping on bipartite lattices. 
The analysis carried out there shows that the positive U Hubbard model away from half-filling {S^ ^ 0) is equivalent 
to a negative U Hubbard model with non-zero J^. Since the above derivation holds independent of the sign of U, it 



4 



follows that [J^] is a smooth, continuous function of (p — (t//2)), as remarked earlier. This implies that the density 
is a continuous function of the chemical potential at any filling, and hence the criterion for phase separation leads 
one to conclude that the one-band Hubbard model does not support phase separation at any filling in one- and two 
dimensions at T = 0. The above conclusions are likely to be unchanged by addition of a two-body term describing 
nearest neighbor interactions at half-filling. The situation away from half-filling requires more care, because the n.n 
two-body term breaks the pseudospin SU(2) invariance away from n=l. However, addition of terms describing next 
nearest neighbor (n.n.n) hopping in the Hubbard model will break particle-hole symmetry at half-filling, and thus 
invalidate the conclusions obtained above. In addition, coupling to lattice vibrations via a local Holstein-like coupling 
might affect our conclusions. Possible extensions of the work presented here to include these interesting cases is an 
open and interesting problem, and is left for the future. 

Alongwith the demonstration of the impossibility of phase separation in the one- band Hubbard model at finite T 
in one- and two dimensions [4], our demonstration of the impossibility of phase separation in the lower dimensional 
one-band Hubbard model has interesting consequences if phase separation is intimately related to the physics giving 
rise to high-Tg. It has been suggested that phase separation (PS) in cuprate superconductors is a relevant issue 
[12,13,14] and may favor superconductive pair formation near the insulator-metal phase boundary. We could say 
that the above derivation implies that it is inappropriate to use the one-band Hubbard model to address the issue 
of phase separation in cuprate SCs. This argument might not apply if the Hubbard Hamiltonian is augmented with 
special terms e.g, n.n.n hopping terms. In this case, particle-hole symmetry is inoperative, and precludes the use of 
formal mathematical arguments leading to EPS. The same argument holds for the extended Hubbard model (with 
nearest neighbor interactions) away from half- filling, where SU{2) symmetry is broken. It should also be remarked 
that the phase separation may probably be of electronic origin, rather than one driven by electron lattice coupling 
[15]. Evidence of charge ordered phases [15 ] near the PS phase could be interpreted as arising from coupling of lattice 
vibrations to a phase separated phase which has its origin in a purely electronic mechanism. 

On the other hand, phase separation arising from a purely electronic mechanism has a natural explanation in models 
with degenerate ground states. To be specific, the three band Hubbard model [16,17,18] (references in [1]) or the 
t — t' — U Hubbard model [19] with/without phonons do show numerical signatures of phase separation. Both models 
break particle-hole symmetry present in the simplest one-band Hubbard model, rendering the above proof invalid. 

Phenomena akin to EPS have also been observed on mesoscopic length scales in the colossal magnetoresistance 
(CMR) manganites [20], and discussed theoretically in terms of effective models. It has been argued that the mod- 
ification of these mesoscale, phase separated regions by external magnetic fields holds the key to the spectacular 
sensitivity of the transport properties of CMR materials to small perturbations. Interestingly, we see that manganites 
require a description in terms of multi-orbital Hubbard models (with ordered/disordered Jahn- Teller distortions). The 
multi-orbital Hamiltonian relevant for describing manganites is not SU (2) invariant in the orbital pseudospin sector, 
but has an effective Z2 symmetry. Thus, we expect a generalised Falicov-Kimball type of model to describe such 
materials [21]. Hence, EPS found in manganites is in conformity with our result here, to the extent that the realistic 
model used for describing their physical behavior is not SU{2) invariant from the outset. Indeed, Freericks et al. [22] 
have rigorously provc^d EPS for the Falicov-Kimball model for any dimension. Rcf. [21] finds a ground state degeneracy 
(within d = 00) resulting from the mapping onto a Falicov-Kimball (electronic version of a binary alloy) model. This 
would imply a jump in ri(/i) as a function of /i, and EPS as a consequence. In fact, such a model exhibits the electronic 
analogue of alloy phase separation, extensively studied in context of alloy physics [23], Finally, this also implies that 
multi-orbital Hubbard models should generically exhibit electronic phase separation: given that insulating, orbital 
ordered (00) ground states in these models spontaneously break the Z^r symmetry of the Hamiltonian {N is the 
number of degenerate orbitals), doping such models may induce EPS close to the correlation driven, first-order Mott 
transitions accompanied by melting of the 00 insulating states. In the CMR case, this is manifested as mesoscopic 
hole-rich (Jahn- Teller (JT) undistorted) clusters in the background of JT-distorted, hole-poor regions. Their size 
increases with modest external magnetic fields because the double-exchange induced increase in carrier kinetic energy 
overcomes short-range 00, and this has indeed been proposed as a mechanism of CMR [20]. 

In conclusion, we have shown that the one-band Hubbard model does not show magnetic LRO in Id, but our 
calculation yields only an upper bound on the sublattice magnetization in 2d. Hence, we conclude that AF LRO is 
not excluded in 2d at T = 0. We have used the fact that the sublattice magnetization is a smooth function of the 
applied magnetic field, along with a particle-hole transformation, to show that the Hubbard model does not support 
phase separation at any filling at T = in one- and two dimensions. Extensions of the derivation to consider the effect 
of special terms like those describing n.n.n hopping and coupling to phonons (where EPS might indeed be realised) 
are being studied and will be reported in future. 
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